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The degradation in bit error rate performance due to imperfect subcarrier tracking J 
by the Demodulation Synchronization Assembly (DSA) is investigated . Results apply to 
any type of digital loop and received signal dynamics. A type four loop causes the least 
amount of loss, because it tracks phase jerk with zero steady-state error . However, when 
I f and f are as large as in the extended Magellan mission, it will be necessary to decrease 
the loop update time in order to minimize the losses, -figures 2 through-#-, illustrate j 



I. Introduction II. Discussion 


Figure 1 depicts, in block diagram form, the overall process 
of demodulating, synchronizing, and decoding a stream of 
binary data. Subcarrier demodulation and symbol synchroni- 
zation are performed in the DSA of the Baseband Assembly 
(BBA). From there, the convolutionally encoded data go to 
the Maximum Likelihood Convolutional Decoder (MCD) for 
decoding. Phase jitter and phase error due to Doppler in both 
the subcarrier and the symbol synchronization loops decrease 
the energy-per-bit to noise spectral density ratio ( E b /N 0 ) at 
the input to the MCD. This decrease in E b jN 0 increases the bit- 
error rate (BER) at the decoder’s output. Given a desired BER, 
the increase in dB of the E b /N 0 necessary to compensate for 
this degradation is denoted as demodulation loss. 

The degradation in BER due to the effects of phase jitter 
and Doppler in the subcarrier tracking loop is estimated. The 
analysis applies to any digital loop with an integrate-and-dump 
circuit. The numerical results, which are applicable to the 
existing and potentially useful BBA’s loops, are summarized in 
Figs. 2 through 8. The subcarrier demodulation losses were 
estimated at a nominal BER of 5 X 10 -3 . 


First, it is necessary to express the variance of the phase 
error as a function of the input signal-to-noise ratio. The 
variance of the open-loop error signal at the output of the 
integrate-and-dump is given by [Ref. 1] 





(i) 


where 


S = the signal power 
K = the number of symbols per update 
K 2 , K 2 , and K $ = gains defined in Ref. 1 
Also, the number of Nyquist samples per symbol is given by 
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where 


B n = the baseband noise-equivalent bandwidth 
r = the symbol rate 

Finally, the thermal noise power is shown as 

< = N oK ( 3 ) 

Assuming that the spectrum of the error process (Fig. 3 of 
Ref. 1) is wide in relation to the loop bandwidth, then the 
variance of the steady-state error signal will be 


(3), (5), (7), and (8) into Eq. (6) and simplifying, the variance 
of the phase error becomes 



With no phase error in the carrier, subcarrier, or symbol 
synchronization loops, the bit error probability for the con- 
volutional code can be expressed by the equation 




Here H(z) is the closed-loop transfer function. Using 
Eq. (82) of Ref. 1 , Eq. (4) can be written as 

a 2 = a 2 2T B, (5) 

e L 

where 



where C x and C 2 are constants which depend on the rate and 
constraint length of the code. The presence of an instanta- 
neous phase error 0(f) in the subcarrier loop degrades E b /N 0 
by the factor 



T = the loop update time 


Assume that 0(f) is of the form 


B l = the one-sided noise -equivalent loop bandwidth (see 
Table 1 of Ref. 1 for various values of B L ) 

Given af w> the variance of the phase error (a^) at update 
instants is obtained from the relation 
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where Gq is the “gain” of the integrate-and-dump device and 
is given by [Ref. 1 ] 
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where 


4> r = the random component of <p(t) and is modeled as 
a Gaussian random variable with zero mean and 
variance given by Eq. (9) 

<j> d (t) = the deterministic component of 0(f) which is of 
the form 

<P d (t) = <t> ss + ^ q t + *f t2 + n f f3 / 3 ( 13 ) 


where 


It can be shown that 


<t> $s = the steady-state phase error at loop update instants 


2 A/L B_ 


( 8 ) 


where E b /N 0 is the ratio of the energy -per-bit to noise spectral 
density. The factor of 2 comes because the rate 1/2 convolu- 
tional code has two symbols per bit. Inserting Eqs. (1), (2), 


0^ = the DCO’s phase rate quantization error 

• •• 

/,/ = the frequency rate and frequency acceleration of the 
received subcarrier which appear because of the 
Doppler effect 

It is assumed that /and higher derivatives are negligibly small. 
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Let f T be the frequency of the transmitted subcarrier, v(t) 
the instantaneous radial velocity of the spacecraft relative to 
the receiving station, and c the speed of light. Then, the 
instantaneous frequency of the received subcarrier will be 



(14) 


where 


p. and z. = the poles and zeros of the loop filter (p.^ 1 1) 

g = normalized computation time 

G = effective loop gain defined by Eq. (9) in Ref. 2 

Assuming that 3a^ « tt/2, the bit-error rate [Eq. (10)] 
together with the degradation factor [Eq. (11)], averaged 
over one update interval, will be 
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Meanwhile, / and / will be the first and second derivatives of 
Eq. (14) 
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which, for v(r) « c, reduces to 




M 

r c 


and 


/ = -/, 


v(0 


r c 


/ 1 “ p [ c 4( 




tt/2 ) 


2a; 


1 — f | exp « (l -^-Vlerlc/ — — 

irJL \ ’ /2 > ( 


d0 r dr 

(18) 

<P d (t) 


2T(d) 


eX P 


/ <^ d (0\2 / \ 
a(l+-^r- erfc / 6 + — — \ 

\ »/2 /J ( (MT.J 




(19) 


where 


d N„ 


b D C *N 0 7T 


The steady-state phase error <t> ss depends on the number of 
integrators present in the subcarrier tracking loop. Using 
Table 2 anci Eq. (30) of Ref. 2, Table 1 is obtained. In this 
table, 
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and <t> d (t) is defined by Eq. (13). 

For the (1/2, 7) convolutional code, C x = 85.7469 and 
C 2 = -5.7230. The second integration over the time variable t 
must be performed numerically. By measuring the horizontal 
distance at a given value of P B , the corresponding loss of the 
subcarrier loop is obtained. 
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III. Numerical Results 

The average bit -error rate ( P B ), expressed by Eq. (19), was 
calculated for several loops of types 2, 3, and 4. The corre- 
sponding loss was estimated at P B = 5 X 10 -3 . 

Curves of subcarrier demodulation loss versus / are shown 
on Figs. 2 through 8. It was assumed that, for the Voyager 
encounter, / would be less than 1.0 mHz/s, and / less than 
1 X 10' 7 Hz/s 2 . For the Magellan mission, it was assumed that 
/ would be less than 25 mHz/s, and / less than 3 X 10 -5 Hz/s 2 . 

The loss for a type 2 loop is quite significant, due to a high 
steady-state error even with a small /. The loss can be reduced 
by decreasing the update time (Fig. 2). 

If a type 3 loop is selected, the loop is degraded only when 
/is significant, which increases the steady-state error. 

In a type 4 loop, the steady-state phase error is insensitive 
to either / or /. However, when / and / are high, as in the 
Magellan mission, (p d (t) can be quite large at the end of the 


loop update instant. This produces a large demodulation loss 
(Fig. 8). In this case, the loss can be reduced only by decreas- 
ing the loop update time. 

In the numerical calculations, it was noticed that the contri- 
bution to the demodulation loss due to thermal noise was 
insignificant for E b /N 0 > 2.0 dB and B L < 1.0 Hz. For this 
reason. Figs. 2 through 8 apply to data rates between 10 and 
500 kilosymbols per second. 


IV. Conclusion 

An equation for the average bit error rate versus E b /N 0 and 
spacecraft dynamics was derived. A type 2 loop is very sensi- 
tive to / and /. A type 4 loop has a zero steady-state error even 
with high values of /. However, with high /and /, the degrada- 
tion due to phase error in the subcarrier loop could be quite 
significant if the loop update time is large. Reduction of loop 
update time will be necessary for the Magellan mission. Fig- 
ures 2 through 8 numerically illustrate the above conclusions. 
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Table 1. Steady-state phase error (at loop update Instants) 


Loop 

Type 

Number of Integrators 
in the Filter 

Steady-State Phase-Error, (t> ss 

f Hz 

/Hz /sec 

/ Hz/sec 2 

1 

0 

V.fTF c 
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- 

2 
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0 
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oo 

3 

2 
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3 
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r BASEBAND ASSEMBLY I 



Fig. 1. Demodulation and decoding process 


DATA RATE = 10- 500 Ksymb/s 
f = 0.0 



L. H. LOOP 

DATA RATE = 10-500 Ksymb/s 
f = 0.0 


Fig. 2. Subcarrier demodulation loss versus f for loop type 2 
at: B l = 0.38, T = 2.65 and 1.0, p 1 = -0.176, p 2 = -0.68, z 1 = 
0.828, z 2 = 0.00, z 3 = 0.0, z 4 = 0.0, G = 0.188, f'= 0.0 


Fig. 3. Subcarrier demodulation loss versus f for loop type 2 
at: B l = 0.19, T = 2.65, p 1 = 0.0, p 2 = 0.0, z 1 = 0.895, z 2 = 
0.1048, z 3 = 0.0, z 4 = 0.0, G = 0.1248 


L. H. LOOP 

DATA RATE = 10-500 Ksymb/s 
f = 0.0 


Fig. 4. Subcarrier demodulation versus f for loop type 2 
at: B l = 0.41, T = 2.65, p 1 = 0.0, p 2 = 0.0, z, = 0.853, z 2 = 
0.1464, z„ = 0.0, z A = 0.0, G = 0.1665 










Fig. 5. Subcarrier demodulation versus f for loop type 3 
at: B l = 0.51, T = 2.65, p 1 = 0.0, p 2 = 0.0, z 1 = 0.217, z 2 = 0.007, 
z 3 , z 4 = 0.883 ± j 0.104, G = 0.1873 
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Fig. 7. Subcarrier demodulation versus f for loop type 3 
at: B l = 0.5/7, 7 = 1.0, p, = -0.173, p 2 = -0.999, z, = 0.960, 
z 2 = 0.960, z 3 = 0.0, z 4 = 0.0, G = 0.15 
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Fig. 6. Subcarrier demodulation versus f for loop type 3 
at: B l = 0.5/7, 7 = 2.65, p, = -0.173, p 2 = -0.999, z, = 0.960, 
z 2 = 0.960, z 3 = 0.0, z 4 = 0.0, G = 0.15 



Fig. 8. Subcarrier demodulation versus f for loop type 4 
at: B t = 0.5/7, 7 = 2.65 and 1.0, p, = -0.173, p 2 = -0.999, 
Z 1 = 0.960, z 2 = 0.930, z 3 = 0.930, z 4 = 0.0, G = 0.15 
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